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Abstract 



Taking as a starting point the results of LDA calculations, which show that in 
MgB2 the phonons have a strong quartic anharmonicity and that the bond- 
stretching electron-phonon interaction (EPI) has both a linear and a large 
quadratic component, we propose a model Hamiltonian which succesfully 
matches a number of experimental evidences. We relate the single critical 
temperature for both superconducting gaps to a phonon-induced inter-band 
coupling whose amplitude increases with temperature. We also obtain phonon 
frequencies and linewidths depending on the band filling, as well as band en- 
ergies and hybridization amplitudes depending on the phonon number. 
PACS:74.20-z, 74.70-b 

I. INTRODUCTION. 

The electronic structure of the 40 K superconductor MgB2 is characterized by the 
presence at the Fermi level of two hybrid bands (cr and tt) of very different character 
[1]. This feature reflects itself in the experimental evidence of two different gaps [2-6], 
which however, in the absence of magnetic fields, have a common critical temperature [7]. 
The observation of a large Boron isotope effect [8] rules out the applicability of theories 
of Coulomb-interaction-driven two-band superconductivity [9], suggesting instead that the 
electron-phonon-interaction(EPl) is the key factor. 

According to the standard theory of the EPI-driven two-band superconductivity [10], a 
single critical temperature for both gaps implies an interaction between the bands contribut- 
ing to the Fermi surface. The microscopic origin of this interaction for the a and tt bands in 
MgB2 is not yet clarified. Impurity scattering can be ruled out [11,12] and, to the best of our 
knowledge, no other precise suggestion has been advanced for the EPI scenario. An estimate 
of the interband coupling strength has been given in Ref. [13], based on the band structure 
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calculation of Ref. [14], yielding a small, but decidedly non-negligible value. The strong 
temperature dependence of the anisotropy of the critical field [13] indicates that the inter- 
band coupling increases with temperature. The present work suggests that such coupling 
might be due to the unusual phononic structure of MgB2. There is a general consensus that 
the dominant electron- phonon interaction is due to a modulation of the inter-site hopping 
amplitudes due to the bond-stretching vibration of the Boron ions [14-19]. By working out 
the corresponding Ehashberg's A , Ref. [15] shows that there is a good agreement with the 
results following the LDA data, and with the experiments. In one-dimensional materials, 
this type of interaction is usually termed the Su-Schrieffer- Heeger (SSH) interaction, and 
we' 11 use this terminology also in the present context for conveniency. 

The unusual phononic features (first revealed by ab-initio calculations of the electron and 
phonon band structures [14,16-18,20]) are the presence of anharmonic contributions (up to 
fourth order in the displacements) to the phononic Hamiltonian, and of both a linear and 
a quadratic term in the SSH interaction. Experimental evidence of anharmonicity comes 
from neutron [16] and Raman scattering [21] data. More specifically, the first-principle 
calculations in Ref. [20] find that the branch has, along the V — A line in the Brillouin 
zone, an energy around 120 meV in the harmonic compound AIB2 and of only 70 meV in 
MgB2, in good quantitative agreement with the neutron scattering data [16]. Additionally, 
Ref. [21] presents first- principle calculations of the evolution of the phonon spectra when Al 
substitutes Mg agreeing with Raman data, which confirm the frequency softening on passing 
from harmonic AIB2 to anharmonic MgB2- One important aspect of such measurements 
is that, as Al substitution changes the occupancy of the bands at the Fermi surface, those 
data show that the phonon frequency and lifetime both depend on the band filling. So, in 
relating those data to the change in phononic properties, one should be able to disentangle 
the effects of anharmonicity from those of band-filling variations. However, one must also 
mention that Raman and infrared data of Ref. [19] for Mgi^xAlxB2 when < x < 0.50 , on 
the contrary, find the frequency of the E2g mode almost insensitive to the Al content. 

Another theoretically predicted effect of the anharmonicity is the reduction of the aver- 
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aged electron-phonon coupling, as expressed by Eliashberg's A. For that, the experimental 
evidence is not so clear. Indeed, the first-principle calculation of the phonon spectra of Ref. 
[18] are in excellent agreement with the experimental data, and they yield at the same time 
that the strength of A is reduced in the anharmonic case. However, if one takes the phonon 
linewidth as roughly expressing the combined intensity of the electron-phonon interactions, 
both Ref. [19] and [21] show that it strongly decreases with x. As Eliashberg's A for a bond- 
stretching interaction [22] depends on structure, band-filling and frequencies, it is difficult 
to precisely identify the cause (s) of the observed effects. 

The ab initio numerical calculations [14,16-18,20] have yielded valuable insights about 
the electronic and phononic structure of MgB2, which we take as the input information for 
the work presented here. Our aim in this paper is to propose a model Hamiltonian which 
represents the physics implied by the results of the first-principle calculations as far as the 
phononic features are concerned. We have no ambition of giving detailed quantitative results. 
However, our model is quantitatively consistent with the numerical results of Refs. [14,16,18]. 
While suggesting a plausible mechanism for the inter-band coupling, and therefore justifying 
[10] the presence of a single critical temperature for both superconducting gaps [2-6], at the 
same time it also qualitatively allows for other experimentally observed features: the increase 
with temperature of the inter-band coupling [13] and the fact that the phonon frequency 
and linewidth both depend on band- filling [19,21]. In particular, we find indications that 
the frequency hardening on Al substitution can not be accounted for by the anharmonic-to- 
harmonic change only. 

Detailed quantitative estimates based on the proposed model are left for future work. 

II. THE ELECTRONIC HAMILTONIAN. 

Our model of the electronic structure of MgB2 by a Hamiltonian has two bands, la- 
belled c and d, which hybridize through an inter-site hopping term. Then, in the real space 
representation, the bare electronic Hamiltonian reads: 
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la la l{j)a l(j)a 

where Hiij) means summing on the z nearest neigbors j of a given site Z, and then 
on In MgB2 one expects that t^.f^f ^ [1,13,15]. The electron-phonon couphng 
parameters in the SSH scenario are derivatives of the hopping amphtudes with respect to 
the inter-site distance. By using — N'^^"^ Y^j^cl^exp (ikRi) we pass to the reciprocal 
space representation, yielding : 



Hei = E + nl. + E + ^4') nt + E ^tt" {cUka + die,,) 

ka ka ka 



(2) 



where t^^ — z ^ X^^^ t'^J exp ik {Ri — Rj) with — c, d. To diagonalize H^i we express 
the bare operators cl, d^ through the hybridized operators a^, according to: 

cL = «L COS ip,. + PI sin (^fe 4a = -«L sin + 4a cos ^j. (3) 

In the following, any operator O expressed through the hybridized operators will be 
denoted by a tilde O. By choosing 



tan(2(^J = 



2ztf 



(4) 



e%-ei + z {tf - tf) 

H^i is brought to diagonal form H^^ =^ H^i — Y,ka {^k''^ka + ^k''^ka) ' particle 



energies in the hybridized bands given by: 



El 



1 



'e" + £^ + z (tf + t 



1 



(5) 
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(6) 



The a and /? bands, at this stage completely decoupled by the transformation of Eq.3, 
represent in our model the a and tt bands of MgB2. 



III. THE PHONONIC HAMILTONIAN WITH ANHARMONIC TERMS. 



Following [16,18] we assume that the purely phononic Hamiltonian H^y^ for MgB2 has 
to include, apart from the usual harmonic term, also a non-negligible quartic contribution. 
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The anharmonicity of the MgB2 phonon modes has been analyzed in [17], showing that 
two modes, degenerate at the F point, have anharmonicities differing in the presence 
(mode labeled a in [17]), or absence (mode b) , of a third-order term. In Ref. [14] a cubic 
term was included in the development of the deformation energy. Its amplitude was found 
to be of the same order as the linear term, i.e. about five times smaller than the quartic 
term amplitude. It gives rise to terms non conserving the phonon numbers, analogously to 
the hnear SSH term. We'll consider only the mode labeled E2g{b) in Ref. [17], which has 
no third-order anharmonicity, both for short, and because the effect of the dropped terms 
is similar, in amplitude and in type, to the one due to the linear SSH term, which we keep. 
A similar assumption of neglecting the E2g{a) mode has been explicitly [23] or implicitly 
[16,18] made in other studies. 

Under such assumptions, Hpj^ takes the form [24]: 

where M is the Boron mass and is the frequency, at the wavevector q along the 
r — A line, of the optical mode of E2g symmetry. The parameter x^^ expresses the strength 
of the quartic term involving the wavevectors ±g and ±p. In MgB2, from Ref. [16], one can 
estimate 7.8 eV-\ 

By quantizing the phonon field according to the usual relations: 



' (8) 



2Mng ' U 1 V 2M V -« U 1 \^ 2Mn^ 

the harmonic part becomes HQ^ {b\bq -\- . When quantizing the quartic term, we 
neglect the terms with different numbers of creation and destruction operators and keep the 
remaining four-operator products only when diagonal. Namely, we approximate bLqb}jb_pbp ~ 
(Sp g + Sp _q^ ^q^-q 1 where 6^6^ = v^. The quartic contribution then reduces to: 

(^) (^) , + \) + *-,) + K) {"i + i-,) « 



qp 



X 



qp 



qp 



X 



qp 



qp \"/\"/ qp \"/\"/ 

The product u^Up is approximated in the MFA fashion, i.e. u^Up pa ^qi^p) 
(Up) (i/g) . Putting together the constant terms, we can rewrite Eq.9 as: 



(9) 



^ (mp 



+ const. 



+ E^^. {blqbl + b_^b^) 
q 

Adding the harmonic contribution and defining 
we obtain the purely phononic Hamiltonian as: 

HpH = E ^^q^q Co + ^0 + E f^^q (^^-.^J + ^-.^ J E ^qp 

q ^ g L P 



(10) 



(11) 



+ const. (12) 



This form can be diagonahzed by a "squeezing " transformation [25] e' 



exp 



T^qTiq (blgbl-b_gb, 



under the condition that 



tanh (2,7]^ 



1/1 



E 



(13) 

^q ' p \ / 

Notice that Eq.l3 yields 77^ < 0. The diagonahzed Hamiltonian e^Hphe~^ can now be 
written as: 



cosh (2r?J + 2 sinh (2r?J (^) E 



b\bq + ^ ] + const. 



(14) 



By substituting 77^ from Eq.13 into Eq.l4, the renormalized frequency uj^ of the harmonic 
Hamiltonian for the squeezed phonons is written explicitly as: 
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1 - tanh^ (2r]^ 



(15) 



where ^qX^ is the phonon frequency entering the quadratic part of the unsqueezed 
phononic Hamiltonian (see. Eq.l2). It is not the true bare frequency, because, from Eq.ll, 
Xq — 1 yields the contribution from the diagonal part of the quartic terms treated in MFA, 
so that ^qXq already contains some effects of anharmonicity, analogous to those taken into 
account, e.g. in Ref. [16]. Thus, Eq.l5 shows that uj^ is increased (hardened) with respect 
to the "harmonic frequency" by > 1., but the squeezing effect, taking account of the 
two-phonon terms terms previously [16] neglected, counteracts the hardening. According to 
Refs. [14,16,18] however, the squeezing effect is not strong enough for an overall softening 
to result. 



IV. THE LINEAR ELECTRON-PHONON INTERACTION. 



The linear part of the SSH electron-phonon interaction is written, in the real space 
representation symmetrized with respect to the site indexes, as 



(16) 



where g^"^ = dt^^ /d (m^ — Ujj \q = —Qji' with ii^u = d, are the coupling constants, and 
the 1/2 factor avoids double counting. 

The Fourier-transformed form of Eq.l6 is written in terms of g^" — (l/z) Qij^ ^wi'^k- 
Aij) (where A^^ is the vector connecting the sites / and j) which we combine in the definition 



of the coupling strength %k-q according to: 



9k-q + C - 9-ik-q) - 9-k) = ^ E 9i[ {sin [{k - q) ■ \ 

(J) 



— sm 



k-A, 



} = iKk-jLq 
(17) 
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Eq.l7 is the simplest possible form of the bond-stretching electron-phonon interaction 
which includes the relevant physics. It is adequate for a qualitative discussion, but it is not 
good enough for a quantitative study. 

Quantization of the phonons according to Eq.8 leads to: 

ep 



When transformed to the hybridized fermion representation H^^ reads: 

ep 

7^ E \^k%-q'^la'^k-q,a + '^k^k-qPkcr(^k-q,a + ^ k^k-q^\a(^k-q,a + ^k,k-qf^ka'^k-q,a {^-q + ^, 

(19) 



Nkqa 



where the effective couplings are defined as: 
^t%-q = lT,k-q COS <Pk COS (Pk-q + TfeJ-g siu (fi^ siu (fi^^^ - -f^k_g siu (y^ + ip,^_g) (20) 

rfffe-g = ^k,k-q sin sin + it^-q COS cos ip,^_g + 7^ji;_g sin ((^^ + ip^_g) (21) 

Tfc J-g = 7fc'fc-9 cos sin - 7^-9 sin (/Pj^ cos (p^-q + l^k-q cos (v'fc + ^k-^ (22) 

rf,!-? = ^k,k-q sin V'ifc cos - 7g_^ cos (^^ sin ^p^_^ + 7^jt_^ cos ((/^j^ + "Pk-^ (23) 

V. THE QUADRATIC ELECTRON-PHONON INTERACTION. 

According to Refs. [14,16,18], the electron-phonon Hamiltonian has to include also a 
quadratic term, which we write in real space in symmetrized form, as: 
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(24) 



where flj^ = d'^t^ /d (ui — Ujj |o = fj'i^, with = c, d. We also develop (ui — = 
+ — UiUj — UjUi- By defining f^" = z~'^ fij^^^''^^^ and introducing the coefficients 
Fk^, = (^/2) {fr + /r - /r, - /p+J , the Fourier transform reads: 

o-(2) _ J_ V [pec f„t „ I „t „ A I ptifi ^ I ^ 



'^q^k-p-q 



+ 



kpqa 



+ ^ E Fkpq (dl^Cp^ + c\^dp^ + (iLpo-C-fco- + ctpcr(i_ 



kpqa 



■ka I "^q^k-p-q 



(25) 



When quantizing the deformations according to Eq.8, we shall take into account only 
the diagonal terms and those which can be diagonalized by squeezing, by enforcing k = p. 
Then UgUj^_p_g reduces to: 



'^q'^k-p-q 



u. 



{^-q^pk = ^pk^l {h-qhl + \h_q + i/g + + l) 



(26) 



Let us stress that our aim is to show that there are some contributions to which 
provide an effective inter-band couphng. We do not claim to be able to treat all the terms 

in i^^p-*: we just want to select the subset of "hot "terms. Under this approximation, by 
using F^^^ = P-k-q > Eq.25 reduces to: 

~ E -^9 Pk^q^\a^kcr + P^d\^d^^ + F^'^ {d\^Cp^ + c\^dp^ 

kqa 



{blgbl + \b_^ + + + l) 



(27) 



Let' s now pass to the hybridized band picture, through the transformation of Eq.3. If 
we define for short the energies: 



Fkq = 



cos^ + Fti sin^ - sin (2(^, 



(28) 
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sin^ cp, + cos' <p, + F^^ sin {2<p,) 



where 



then Eq.27 can be rewritten compactly as: 



fcgcr 



(29) 
(30) 



(31) 



{hl^hl + + i/^ + + 1) 



(32) 



VI. THE ELECTRON-PHONON HAMILTONIAN IN THE SQUEEZED 

PHONON REPRESENTATION. 

Let us now introduce the squeezed phonon representation also for H'^ + H'^ . By using 
the relation {h'^_q + — e^i {h'^_q + the linear coupling term becomes: 



kqa 



^k%-q(Aa^k-q,a + ^k'!k-qPlaPk-q,a + ^ k^k-q'^kaPk-q,^ + ^k,k-qPla'^k-q,a (^-? + ^. 

(33) 



Then the hnear coupling has a reduce amplitude, as r/^ < (see Eq.l3), consistently with 
the numerical analysis of Ref. [18]. 
For the quadratic part H^^^ we get: 



kq 
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{bl,bl + \h_^ + V^ + + l) 



(34) 



In Eq.34 we can decouple the electronic and phononic terms in MFA. Indeed, assuming 



{h\h\ + \h_^ + l/^ + + 1) 2{v^ + 1 

and defining 



(35) 



(36) 



yields: 



+ ^ 6^'''^$^ + + i/g + + l) + const. 



(37) 



By writing Eg ('^g + |) = | Eg ('^c 
in the hybridized basis we get: 



1 ' 



1^ and reordering the Hamiltonian 



e'Ee-' = E Elnl^ + J] nL + E H + ^) + ^"^i^"" + (38) 



Inserting e^H^e ^ from Eq.37 into Eq.38 and reordering yields: 



e^i^e-^ = Tifr+^ef +^p?. + (^'H^J(^~' + const. (39) 



where the diagonal electronic Hamiltonian is: 



kcr . 



<. + E 

fecr - 



^f + ^E<e^^^ (2H) + i 



n' 



(40) 



and describes a phonon-depending renormalization of the band energies. The hybrid 
electronic Hamiltonian is: 



(41) 
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It represents a phonon-depending band hybridization term. This is, we beheve, the term 
responsible for the coupling between the bands in MgB2 which results in a single critical 
temperature for both gaps. We would like to point out that, in the limit of small inter-bare- 
band hopping, i.e. tij/tfj 0, and at zero temperature, where (i/g) can be neglected, one 
finds: 

J^^, - SL^e^-^F.t = 2Lje2"'. ^ {^os (^'A,,.) [l - cos (gA,,.)] } (42) 

h ' Ij^ <J) 

Apart from geometric factors, this amplitude depends only on the intensity of the squeez- 
ing (through e^^i) and on the strength of the quadratic inter-band SSH electron-phonon 
couphng fif. As fif is a second derivative of tfj, it can be non-neghgible even if tfj itself 
is very small. Different evaluations of flf [1] all agree that in MgB2 it has an appreciable 
value. More specifically, an effective two-band model derived from first-principle calculations 
[14] yields X^J = 0.21 and X^jKa = 0.15 [26]. 

At nonzero temperatures, the hybridization amplitude gets an additional contribution 
proportional to 2(1/^), hence it increases with temperature, consistently with the findings of 
Ref. [13] 

The purely phononic term 

'Hph ^ E (^H + e'"^*.) + + l) + E ^""'^ {bl,bl + b^b_^) (43) 

is diagonalized by a second squeezing transformation = exp — Y,q '^q {^-q^l ~ ^-q^c 
with the value of 1?^ set by: 



tanh (2d^ = ^ - 1 (44) 



Notice that the sign of 1?^ is opposite to the sign of Also, as the relevant phonons 
are optical ones, hoUg never vanishes in the Brillouin zone, then tanh (^-dgj ^ —1 , and -d^ is 
always well defined. Due to the presence of $g in Eq.44, the parameter 1?^ depends on the 
band- filling factors (n'^i^^) and on the band hybridization (q;|.^/3^^ -|- Pl^ai^^). 

The diagonalized free-phonon Hamiltonian reads therefore: 
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e^npf,e~^ = J2 [^^q cosh (2i?J + 2e^''«e^''«$J + const. = ^ /izu^i/^ + const. (45) 



We have obtained a band-fiUing and hybridization-depending renormahzation of the 
phonon frequencies: 



hvjg = ;iu;gCOsh (2^?J + 2e2''<'e^''«$g 



(46) 



The final Hamiltonian therefore reads: 



kcr kqa 

+ ^ /iG^gi^g + H^^ e~'^ + const. 



(47) 



where 



4 -Ei + ^Y. Fiy^^ (2(^,) + l) (C = «, /3) 



(48) 



This has the shape of a standard (i.e. harmonic and hnear) SSH Hamiltonian for two 
hybridizing bands. However, also the linear SSH term has acquired a band-filling and hy- 
bridization dependence, because it now reads: 

e^'''^^'^ ^t%-q^ka^k-q,a + ^k'^k-qPkaPk-q,a + ^ k'^k-q^kaPk-q,a + ^k,k-qPka'^k-q,a (^-( 



kqa 



(49) 



Therefore also the phonon linewidths, due to the SSH interaction, will depend , through 
exp(i?J, on (riki^^) and (4^/3^^ + ^^a^^) . 

The weakening of the linear electron-phonon interaction is expressed by the coefficient 
exp(i?^) exp(?7g) . The value of exp(?7g) is set by the diagonalization condition of the anhar- 
monic phonon Hamiltonian, Eq.l3 according to the identity: 



1 + tanh 2r]^ 



\1- tanh (27]gj 



(50) 
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Therefore the squeezing effect related to the anharmonicity of the phonons also reduces 
the electron-phonon interactions. 

To conclude, let us check if the link that our model establishes between the renormal- 
ization of the harmonic frequency from fi^ to cUg and the reduction of the electron-phonon 
coupling strength is consistent with the estimates of those quantities as given, e.g., in Refs. 
[14,16,18]. The ratio /w^is evaluated as 85%. [14], 75% [16] and 80% [18]. We assume that 
the value of exp(i?g) (Eq.44) , which we can not estimateat this stage, is not far from unity. 
By taking from Ref. [16] x^g ~ 7.8 eV~^, hQg — 60 meV, and assuming a dispersionless mode 
in Eq.ll we obtain, at zero temperature, where we can neglect {Ug), that ~ 1.12 . Then 
tanh (2?7g) ~ -0.105 yielding, from Eq.l5, Qq/cOg ~ 0.85 and, from Eq.50, exp (277J = .90. 
In Ehashberg's theory [22] A ~ {\g'^\) / {uj'^) , where (...) are suitable averages. In our model 
then Xanhar/ ^har ~ c^'^'j / 1 — tanh^ (^^9) ~ 0.65. which agrees with the estimate [18] 
of a 30% weakening of Xanhar- Prom the fact that tanh (277^^ ~ —0.105 we also conclude 
that the large phonon softening found on passing from AIB2 to MgB2 [16,20,21] is probably 
not due only to the harmonic-to-anharmonic phonon change, because the squeezing-induced 
softening effect in Eq.l5 is not strong enough. The filling-dependent effect of exp('j?q) should 
also be taken into account quantitatively before drawing more reliable conclusions on this 
point. 



VII. CONCLUSIONS. 

We have obtained a model Hamiltonian which should contain the essential physics of 
MgB2. 

Our starting point was a two-band Hamiltonian with anharmonic phonons and both 
linear and quadratic electron-phonon interactions of the bond-stretching type, as dictated 
by the results of LDA calculations [14,16-18,20]. The final Hamiltonian followed by ap- 
plying a sequence of unitary transformations to both the electronic and phononic terms. 
In particular, we have been able to go beyond the Hartree-Fock approximation in treating 
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the anharmonic effects. We have thus obtained an effective Hamiltonian (Eqs.47 and 49) 
of a very simple structure. The electronic part has two bands with a phonon-depending 
hybridization, generated by the quadratic electron-phonon interaction, which increases with 
temperature as observed [13]. The phononic part has an effective harmonic free-phonon 
term with a frequency which depends from the band filling factors. Finally, the effective 
electron-phonon interaction is reduced to a linear one, but with an amplitude also depend- 
ing from the band filling factors, which would result in a similar dependence of the phonon 
hnewidths [19,21]. One could ask if it might have been possible to obtain the same results 
starting from a simpler Hamiltonian, as, for instance, the one proposed in Ref. [23]. We do 
not think so. Indeed, all our results depend basically from taking into account the quartic 
anharmonicity and the second order EPI, as can be checked by considering the limits for 
?7g — > and — >^ of Eqs.46,47 and 49. Namely, if both rjg, i!}g — > 0, from Eq.49 one finds no 
reduction of the effective EPI. If there is no second-order EPI, then = follows, implying 
•dq — > 0, SO that there is no filling dependence of zUg (Eq.46), no phonon-number effect on 
the hybridization (Eq.47) and no filling dependence of the phonon linewidths (Eq.49). 

We have also shown that the numerical results of Ref. [18] about the phononic features 
of the material, namely the renormalization of the effective harmonic frequency oUg and 
the reduction of the Eliashberg's A, can be consistently interpreted as due to the phonons 
accomodating themselves in a "squeezed "state. On the other side, squeezing effects alone 
are not strong enough to account for the large phonon softening on passing from AIB2 
to MgB2 [16,20,21]. 

Acknowledgement 1 We thank J. Spalek, M. A. Gusmdo and J. R. Iglesias for stimulating 
discussions. 
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